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We develop a diagrammatic approach with local and nonlocal self-energy diagrams, constructed 
from the local irreducible vertex. This approach includes the local correlations of dynamical mean 
field theory and long-range correlations beyond. It allows for example to describe (para-)magnons 
and weak localization effects — in strongly correlated systems. As a first application, we study the 
interplay between nonlocal antiferromagnetic correlations and the strong local correlations emerging 
in the vicinity of a Mott-Hubbard transition. 
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I. INTRODUCTION 

Strongly correlated electron systems represent both an 
opportunity and a challenge for modern physics: An 
opportunity, since fascinating phenomena occur such as 
high-temperature superconductivity in cuprates, "colos- 
sal" magnetoresistance in manganites, and quantum crit- 
ical behavior in heavy Fcrmion compounds. But at the 
same time they are a challenge, since the very same corre- 
lations which are responsible for these phenomena make 
a theoretical understanding and hence an experimental 
optimization of these effects particularly difficult. 

One of the key issues which arises due to strong elec- 
tronic correlations and which cannot be described by per- 
turbation theory is the Mott-Hubbard metal-insulator 
transition.^ In this respect, dynamical mean field the- 
ory (DMFT)^-"^ was a big step forward to a more thor- 
ough understanding of this transition. DMFT becomes 
exact in the limit of high spatial dimensions {d — > oo) 
and already accounts for a large (local) part of elec- 
tronic correlations — the part which provides for the radi- 
cal changes upon going from a metal to a Mott-insulator. 
Real physical systems are however one-, two-, or three- 
dimensional. Hence, nonlocal correlations, which are ne- 
glected in DMFT, may be of importance. Corrections of 
order 1/d have been considered in Ref. 0, resulting in a 
two-impurity problem, and account for short-range cor- 
relations. There has also been recent progress to go be- 
yond DMFT through cluster extensions)^ which include 
correlations within the cluster. These correlations are 
also necessarily short-range in nature due to numerical 
limitations of the cluster sizei^ 

Often, however, long-range correlations are of vital 
importance. They are responsible for a rich variety 
of phenomena, ranging from magnons and screening of 
the Coulomb interaction to quantum criticality. Long- 
range correlations are also generally pivotal in the vicin- 
ity of phase transitions. The existing theories describing 
long-range (e.g., magnetic) correlations such as the fluc- 
tuation exchange approximation^, the two-particle self- 
consistent approximation*, and the functional renormal- 
ization group^ are restricted to the weak-coupling regime. 
For strongly correlated systems, e.g., in the vicinity of 



a Mott-Hubbard transition, an extension of DMFT by 
nonlocal (particularly long-range) correlations is hence 
needed. 

For static mean-field theories, such corrections have 
been studied since decades, e.g., for localized^? and itin- 
erant magnetSfiiii^ for disordered systems such nonlocal 
effects have also been considered ^13^ But there have been 
only a very few attempts so far to include long-range 
correlations beyond dynamical mean field theory: The 
DMFT self-energy was supplemented by an "external" 
k-dependent self-energy which describes spin fluctuations 
in the spin-fermion niodel^'* or which stems from the self- 
consistent renormalization theory, and one might also 
subsume the GW-I-DMFT approach^^ here. Let us also 
note the extended DMFT (E-DMFT^"^), which consid- 
ers the effect of non-local interactions on the purely local 
self-energy. 

In this paper, we aim at a systematic diagrammatic 
extension of DMFT by long-range correlations and at 
an investigation of their effect on the non-local self- 
energy. Diagrammatically, DMFT corresponds to all 
topologically distinct, but local Feynman diagrams for 
the self-energy. On the next level we assume the locality 
of the fully irreducible two-particle vertex, and consider 
all (local or nonlocal) self-energy diagrams which can be 
constructed from this vertex. One might generalize this 
approach, requiring locality of the fully irreducible n- 
particle vertex. Then, one has DMFT for n = 1, the 
dynamical vertex approximation (DFA) for n = 2, and 
the exact solution for ti ^ oo. We think however that 
the one- and two-particle levels [n = 1,2) are the most 
relevant approximations. If one is interested in a specific 
physical problem, a restriction of DFA to certain lad- 
der diagrams is reasonable. In the particle-hole channels 
the ladder diagrams yield (para-)magnons^^ and RPA 
screening; in the particle-particle channel the cooperon 
diagrams are responsible for attractive pairing interac- 
tions and weak localization effects. DFA includes such 
ladder diagrams beyond DMFT, but with the local vertex 
instead of the bare interaction so that strong correlations 
are accounted for. 

In this paper we introduce DFA and apply it for 
studying long-range antiferromagnetic fluctuations in the 
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three-dimensional Hubbard model. The interplay of 
these nonlocal spin fluctuations with the local DMFT 
fluctuations is surprising: Close to the metal-insulator 
transition the nonlocal fluctuations strongly suppress the 
spectral function, in contrast to the weak-coupling expec- 
tation for three dimensions. 

The plan of the paper is the following: In Sec. II, we 
introduce the DP A. Specifically, we derive the full DP A 
scheme based on the parquet equations in Sec. IIA, and 
a simplified version based on a ladder subset of diagrams 
in Sec. TIB. The latter yields the most important dia- 
grams for the specific problem considered in this paper, 
i.e., the paramagnon fluctuations in the proximity of the 
AF transition. We compare this approach with the \/d 
expansion in Sec. IIC. The details concerning the calcu- 
lation of the local vertex within an exact diagonalization 
impurity solver are reported in Sec. Ill (and in the Ap- 
pendix). Results for the local irreducible vertex and the 
DFA self-energy and spectral functions are presented in 
Sees. IVA and IVB respectively. Finally we give a con- 
clusion and discuss the potential of our new method in 
Sec. V. 



II. DYNAMICAL VERTEX APPROXIMATION 

Starting point of our considerations is the Hubbard 
model on a cubic lattice 



H 



(1) 



where t denotes the hopping amplitude between nearest- 
neighbors, U the Coulomb interaction, c\^{cia) creates 
(annihilates) an electron with spin a on site «, riic — 
c\^Cic,. In the following, we restrict ourselves to the para- 
magnetic phase with n electrons/site and temperature T. 

Let us suppose we know the two-particle vertex 
-'^kk'q^^- Then, we can calculate the (nonlocal) self- 
energy through the exact relation (following from the 
equation of motion; see Fig. la and, e.g., Refs. 20,21) 

Sk.i/ ~ ^ E -'^kk'q^^^k'+q.l^'+wGj^V' Gk+q,y+w(2) 

u' UJ 
k'q 

where Gk,iy = (*t'n~ek+/^~Sk,i/)~^ is the nonlocal Green 
function, ek = — 2i ^^^^ ^ ^ cos Aiq, the bare electronic 
dispersion, and the electronic chemical potential. Gen- 
erally, Fj^J^/q^''' can be expressed diagrammatically, e.g., 
by taking the fully two-, three- and more particle irre- 
ducible local vertices as building blocks and connecting 
these blocks by local and nonlocal Green functions. 



A. Parquet equations 

In the DFA, we restrict ourselves to the local fully- 
irreducible two-particle vertices F|?5. From these build- 



f k'+q,J- A 


f k'+q,i \ 








k"+q,CT 


AW 



kt 



k+q,t 



kt kt k+q,t 



kt 




FIG. 1: a) From the reducible vertex we directly obtain the 
self-energy, b) and c) The two particle-hole channels con- 
tributing to the self-energy (longitudinal and transversal) in 
the ladder approximation. Instead of the bare interaction, 
ladder diagrams are constructed from the irreducible local 
vertices (crosshatched) in DFA. 

ing blocks, the reducible vertices Fj^J^,'^^''' can be ob- 
tained through the self-consistent solution of the parquet 
equations] -^^1^° Representing this vertex as a sum of con- 
tributions of different channels, one hasi^iSSiiSi 
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Here, f;;;;';^^^ = contribution of 

the three channels can be written in the following form 
(see Fig. [2]) 
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(5) 
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^kk'q 



(FTT+FTi)*G*G*(FTii^^ + G + Z) 

- FTUG*G*(f1jii„^ + G + Z) (7) 

where * stands for multiplication and summation over 
the momenta and frequencies given in Fig. [21 

For determining the fully irreducible local vertex F}?'^, 
which is constructed from purely local Feynman diagrams 
only, we resort to the Anderson impurity model (AIM)-^». 
In fact, the AIM has only one interacting site, so it yields 
the same local diagrams -and hence the same irreducible 
vertex- provided that the local Green function is identi- 
cal. Hence a practical way to obtain F|?'^ is through the 
(e.g., numerical) solution of the AIM. Starting point can 
be the (local) spin and charge susceptibility of the AIM 



As(c),loc Aloe (-) '^loc 
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from which we can obtain the full (reducible) local ver- 
tex Ts(c).ioc and the irreducible local vertices in the spin- 
, charge- and particle-particle channel (Fg^jr, F^ir and 
Tpp,!!-, respectively) via the standard relations: 

I^I/'w A , -I- ■^/'^ Y.V}j'uj v' 

As(c),loc — AOwJoc"!^!'' + XOu,\oc^ s(c),locAOa;Joc 

(10) 



-'-s(c),loc — 11^ s(c), ir7 AOw,loc'Viy 



pp,loc 

where 



[(l^p'p.ir )u.u ^t)i>+vWoc^'^i>\o=u+uji'^'^) 



xL.ioc = -TGioc(i'')Gioc(t^ + y'), 

nL.loc=rGloc(^^')Gloc(c^~^^')- 



(12) 
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FIG. 2: Graphical notation 
of the parquet Eqs. ([3]), (O, 
P, and ©. C, Z, and 

Z denote the contributions of 
the particle-particle, particle- 
hole and interaction chan- 
nels to the non-local ver- 
tex (^^e mo- 

mentum and frequency k2 , V2 
are determined by conserva- 
tion laws). 



From these vertices we can in turn calculate the fully 
irreducible vertex as 
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energy, which is not relevant near magnetic instabilities, 
has been neglected here. 

The quantities on the right hand side of Eq. (|14p are 
calculated from the local vertex T^^'^^^j^, irreducible in the 
spin (charge) channel, via 



B. Ladder approximation 

In this paper, we are particularly interested in param- 
agnon contributions affecting the self energy in the vicin- 
ity of the antiferromagnetic phase. Hence, as discussed in 
the introduction and with some justification in l/d (see 
Sec. II C), we restrict ourselves to the ladder subset of 
the parquet diagrams in the two particle-hole channels 
shown in Fig. lb and Ic. These diagrams can be derived 
from the general parquet set of diagrams of Fig. 2, sup- 
posing the locality of C, Z, Z in the r.h.s. of Eqs. ©-(0). 
Expressing the ladder diagrams through the vertices in 
the spin (s) and charge (c) channels ^'^Ji^J^^, which de- 
pend on the momentum transferred q only, the sum of 
the vertices of Fig. lb and Ic is obtained as 



^ s(c),q 



ic), 



(15) 



where Xo'^cu = -TJ2i^Gk,^'Gk+qy+u, with Gk,,. = [w - 
eic + M - Eioc(i^)]"\ Sloe being the local (DMFT) self- 
energy. Note that contrary to the full parquet solution in 
Sec. Ill Al the self-energy of the internal Green functions 
is considered purely local in accordance with the assump- 
tion of the locality of the vertex The results of 
this non-self-consistent approach are expected to be close 
to those of the self-consistent one, due to the cancella- 
tions between (self-consistent) non-local self-energy and 
corresponding corrections to the vertex ^s{c) in cf. Ref. 

n. 

Substituting Eq. into Eq. we obtain after a 
shift of the momenta and frequencies 
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(14) 



Here, the first two terms of Eq. describe the lon- 
gitudinal and transverse paramagnons in Fig. lb and 
Ic, respectively, and the last term subtracts the double- 
counted local contribution. Note that the nonlocal con- 
tribution of the particle-particle channel to the self- 
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Eq. (|T6l) reduces to the DMFT self-energy if the nonlocal 
quantities are replaced by local ones. But beyond that, 
it describes the nonlocal ladder diagrams of Fig. lb and 
Ic. 
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III. CALCULATION OF THE LOCAL 
FOUR-POINT VERTEX 
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FIG. 3: Diagrams for the self-energy in terms of the vertex 
to order l/d°, [diagram a), i.e., DMFT] and [dia grams 

b)-f)]. The ladder diagram g) show how diagrams b)-e) can 
be constructed from a fully irreducible local two-particle ver- 
tex. In contrast, for diagram f) a local three-particle vertex 
is needed, see h). Hence the contribution f) is the only one 
not included in the DFA. 



C. Comparison to the 1/d expansion 



Let us compare the result Eq. (fTB|) to that of the l/d ex- 
pansion. While the DMFT self-energy contains the local 
Green functions and vertex only (Fig. 3a), the leading 
nonlocal corrections to the self-energy are proportional 
to 0(l/(i^/^) and the consideration of the diagrams con- 
taining two different sites is sufficient at this order 4 The 
possible types of these diagrams for the self-energy are 
shown in Fig. [5}D-[5f. The first type of corrections (Fig. 
inb) involve the non-local Green functions only. The sec- 
ond type of corrections (Figs. [3l:-[3^) contains the non- 
local vertices with two legs at % and two at j sites. To 
leading (zeroth) order in l/d these vertices can be ex- 
pressed as a ladder of local vertices connected by the 
non-local Green functions, as in Fig. [3^. The contri- 
butions of Figs. [3Jd-[3^ are hence included in the DFA 
with the ladder approximation, i.e., in Eq. p6p . The 
last type of Xjd^l"^ corrections to the self-energy (Fig. 
|3f) involve the three-particle local vertex (Fig. [3h) (and 
are of order [U/t)^). According to the classification of 
the introduction, these corrections should be taken into 
account on the next level of approximation beyond DFA 
and are the only corrections neglected in DFA. 

Therefore, the DFA reproduces correctly the leading \/d 
correction to the self-energy with the three-particle local 
vertex neglected. 



The calculation of the self-energy requires the 
knowledge of the local vertex, either fully irreducible — for 
the general scheme — , or irreducible in the spin (charge) 
channel — for the ladder diagrams of Fig. lb and Ic [Eq. 
([TS])]. As already noted, this local vertex can be ob- 
tained numerically from the Anderson impurity model. 
For obtaining the four-point vertex ^^gi^^^^^, we need 
to calculate the AIM susceptibility for three Matsubara 
frequencies^^ 



UU UJ(T(T 



Xlo 



l/T 



dTi dT2 dT3 e-'^'" e*^^(''+'^) Q-^r3{^'+^) 



{Tr cl{ri)c.a{r2)){Trcl,{r,)c.a'm] . (17) 



where v, v' and u are the two fermionic and the bosonic 
(transferred) Matsubara frequency, respectively; {T^- ■ ■ ■) 
indicates the thermal expectation value of the time- 
ordered operators and the last term represents the non- 
connected contributions. With a spin (anti-) symmetriza- 
tion ([8]) we obtain the corresponding charge and spin sus- 
ceptibilities Xs('c)^ioc- From these, we can either determine 
the fully irreducible local vertex and through the parquet 
equations the reducible vertex and the self energy along 
the lines of Eqs. ([3]), and ^\ or we can directly 

calculate the particle-hole ladders along the lines of Eqs. 
(ITU) and We implement the latter by (i) solving 

the DMFT equations using exact diagonalization (ED), 
(ii) calculating via Eq. pT|) the local vertices, and (iii) 
constructing from these through Eq. the DFA self- 
energy. In principle, this fc— dependent self-energy yields 
new Green functions and a new vertex. However, because 
of the ladder approximation, we do not perform such a 
self-consistent calculation here. 

Within ED, the calculation of xioc [Eq. p7| ] is straight- 
forwardly (albeit lengthy) performed by resorting to its 
Lehmann representation, whose explicit expression is re- 
ported in the Appendix. The Lehmann representation of 
Xioc requires four summations over all the Hilbert states 
of the discretized AIM (compared with only two sum- 
mations for evaluating the local Green function). This 
is the higher computational cost of DFA compared with 
DMFT. By performing a parallel computation of Xloc^ we 
were able to calculate AIMs with Ns — 5 sites and evalu- 
ate xioc for the lowest N,nax — 20 (or, in some cases, 25) 
Matsubara frequencies. This turned out to be sufficient 
for getting a stable analytic continuation of Eq. (|16p . us- 
ing the Fade algorithm. For the momentum summation 
of Eq. pB|) we have used Nk = 96 points for each direc- 
tions. 
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V/D 

FIG. 4; (Color online) Dependence of the local vertex 
r^"^'"^ " '^ii incoming Fermionic frequency v' , for the 
three different values of U and T indicated as symbols in 
the left inset, which shows the DMFT phase diagram with 
paramagnetic metallic (PM), insulating (PI), and antiferro- 
magnetic (AF) phase. Right inset: Same as main panel but 
at fixed U — ID and at different uu's (circles: uj = 0; triangles: 
uj = 6'kT; squares: uj = 12nT). 

IV. RESULTS 
A. Local vertex close to a Mott-Hubbard transition 

First, we discuss our results for the local irreducible 
vertex. Since we are mainly interested in long-range an- 
tiferromagnetic fluctuations, we consider temperatures 
slightly above the DMFT Neel temperature (see the inset 
of Fig. d]). In this case the largest contribution to the non- 
local part of the DFA self-energy stems from the terms of 
Eq. pH]) proportional to the local vertex irreducible in the 
spin channel, particularly with zero bosonic frequency. 

In Fig. m we show the r^'(j.'"~° as a function of u' (with 
ly fixed to its lowest value ttT) at half-filling for three 
different values of the Hubbard interaction {U — ID, U = 
1.5D and U = 2_D, where D = 2^/6 1 is twice the variance 
of the non-interacting density of states) and temperature 
(T = 0.067D, 0.0891) and O.ID). This choice allows us 
to highlight the remarkable differences occurring when 
moving from the metallic regime to the crossover region 
of the DMFT phase-diagram of the Hubbard model. In 
particular, we note that correctly approaches the 

corresponding value of the bare interaction U for large 
v' . On the other hand, at small we observe a radically 
different behavior of the local vertex depending on the U 
value: At U — ID, Ts.ir displays a smooth minimum in 
the region of small i^' , while at U — 1.5D (and even more 
at [/ = 2D) a very pronounced maximum of Fs appears 
at I/' — V. While the behavior of Vs,ir at small U can be 
easily interpreted as the screening of the bare interaction, 
typical of the metallic phase, the huge maximum of T s^ir 
at larger U stems from particle-hole fluctuations in the 



vicinity of the metal-insulator transition, as discussed in 
Ref. [Hi 



B. DFA self-energy and spectral function 

The striking behavior of F^ has consequences for the 
DFA self-energy on the real axis and spectral function, 
presented in Fig. [5] for the same three different /7's and 
T's of Fig. H (left inset). 

At [/ = ID, i.e., in the metallic regime of the phase di- 
agram, DMFT shows a quasiparticle peak which is only 
weakly damped for k-vectors on the Fermi surface. In 
DFA, the quasiparticle scatters at nonlocal antiferro- 
magnetic fluctuations, resulting in a broadening of the 
quasiparticle with a now significant damping given by 
ImEk(O). For this U value, the distinct features of Ts^ir 
are not yet particularly pronounced in Fig. [H It is the 
strongly enhanced antiferromagnetic susceptibility close 
to the Neel temperature which leads to this damping 
in Eq. Stronger damping effects can be observed, 

studying the Hubbard model iw d — 2, e.g. by means of 
the cluster extensions of DMFT—, which predict a pseu- 
dogap opening at low temperatures. Actually, a stronger 
damping in d = 2 than in the three dimensional case 
considered here is expected from weak coupling perturba- 
tion theory. Unfortunately, to our knowledge, no cluster- 
DMFT calculation has been performed for the case of 
d = i, since it poses severe constraints on the cluster 
size. 

At [/ = 2D, the stronger electronic correlation re- 
flects in more pronounced changes of the spectral func- 
tion. Now, it is the huge Vs,ir of Fig. 5] which strongly 
suppresses the spectral weight at the Fermi level. This 
weight is transferred to the Hubbard subbands, which 
get some additional structure as an effect of magnetic 
fluctuations. The nonlocal fluctuations result in a much 
more insulating solution, albeit the Green function is still 
non-locaP^ . 

We emphasize that the mechanism of spectral weight 
suppression at U = 2D is very different from that in 
the weak-coupling regime (Ref. [l^ and our U = ID re- 
sults) , where long-range magnetic fluctuations in the im- 
mediate vicinity of the magnetic phase transition play 
the key role. In contrast at U = 2D, already relatively 
short-range spin fluctuations are important because of 
the strong correlations reflected in the enhanced Tg^ir. 
The spectral weight suppression is therefore also quite 
robust upon increasing T, i.e., upon going further away 
from the antiferromagnetic transition (not shown). 

At U — 1.5Z), we have something in between the two 
cases discussed above: The vertex is already enhanced, 
but long-range antiferromagnetic fluctuations are still es- 
sential. This leads to a suppression of the quasiparticle 
weight and structured Hubbard bands. Altogether this 
shows that scattering at nonlocal fluctuations close to a 
strongly correlated antiferromagnet is very different from 
that in the vicinity of a weakly correlated Slater antifer- 
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FIG. 5: (Color online) Self 
energy (left: real, middle: 
imaginary part) and spec- 
tral function (right) at k = 
(V2,V2,7r/2), U^ID (top), 
1.5D (central), and 2D (bot- 
tom) and the same T's of Fig. 
H Compared to DMFT (light 
blue line) quasiparticles are 
damped through scattering at 
nonlocal spin fluctuations in 
DFA (dark red line), and the 
system is more insulating. At 
U = 1.5D and 2D, these ef- 
fects are drastically enhanced 
because of strong local corre- 
lations reflected in Fs ir 



romagnet. 



V. CONCLUSIONS 

We developed the dynamical vertex approximation 
(DFA) based on the assumption of the locality of the 
irreducible vertex. For the half-filled three-dimensional 
Hubbard model, we found that the local vertex (irre- 
ducible in the particle-hole spin channel) strongly de- 
pends on all three frequencies; it is hugely enhanced at 
some particular frequencies for large U. These strong 
local correlations entail similarly strong nonlocal fluctu- 
ations. The scattering of the quasiparticles at these non- 
local spin fluctuations, in turn, drastically reduces their 
life times; spectral weight is transferred to the Hubbard 
bands which develop some additional structure. These 
nonlocal effects of strong electronic correlations are very 
different from those at weak coupling. It is a strength of 
DFA to reveal them. 

Including long-range correlations, DFA opens the door 
to study a wide variety of physical phenomena, previ- 
ously described only for weakly correlated systems, such 
as magnons in strongly correlated (anti-)ferromagnet9^, 
the interplay of weak localization effects and strong elec- 
tron interactions, and vertex corrections to the RPA 
screening. A self-consistent realization of the approach 
might also allow us to study how nonlocal fluctuations 
suppress magnetic long-range ordering, whether anti- 
ferromagnetic fluctuations in the vicinity of the metal- 
insulator transition result in unconventional supercon- 
ductivity, and how physical quantities change in the 
vicinity of ferro- and antiferromagnetic quantum critical 
points. 
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APPENDIX: LEHMANN REPRESENTATION 
OF THE LOCAL SUSCEPTIBILITY 

In this Appendix we report the explicit expression 
of the Lchmann representation for the local (and spin- 
dependent) susceptibility, which is necessary for perform- 
ing the calculation of the basic "brick" of the DFA (i.e., 
the local four point vertex) within the ED algorithm (see 
Sec. HI). 

We start with the evaluation of the T-ordered product 
appearing in the definition of the (spin-dependent) local 
susceptibility in Eq. (fT7|) : 

Jo 
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of the equation, which are already ordered in terms of 

decreasing times (with no further sign change). The six 

different contributions Xiocj Xioc^ ''' a-ppearing in the 

last line of Eq. (jA.ip reflect the six different ways of 

arranging the order of the three Matsubara times ti, T2 

-tv'r '^^ ™ integral of Eq. xJoi*, xkH, ■ ■ ■ 

X g iv\Tx T2} gjcj(T2 T3j g 11/ T3 ^ . . . ^ gg^j^ |-,g expllcltly cxprcsscd in a very convenient way for 

ED scheme, i.e., in terms of the eigenenergies Eiq and 

= 123 , 132 , 312 , 213 , 231 , 321n /a 1 ^^t^^^ elements {N\c^J^\M) = (c^^Va/, of the as- 

2^ \X\oc + Xioc + Xioc + Xioc + Xioc + Xioc ) l^-ij sociated AIM, through the standard Lehmann represen- 
tation. The evaluation of the six time integrals in Eq. 
where /3 = 1/T, Z is the partition function and with (|XT|l is straightforward- albeit lengthy, and yields the 
( • • • ) we indicate the thermal and quantum average of following results: 
the four fermionic operators in the r.h.s. of the first line 
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M 



1 



-PE^ 



+ e 



-f)Es 



i{v + Lo) + Em — El \ ilo + En — El iv — En + E 



M 



Xioc 



E 

N,M,L,S 



i(l^' + LU) -E 



M 



-(JEt. 



X {c\)NM{Ca)ML{c\,)Ls{Co')sN (A. 2) 
q-PEl _ g-/3Sjv 



-pEs 



1 



Xioc 



i{v - v') + Em - Es + Em - En iv' + Es - En 
1 1 



i{iy + Lu) + El - Es \ iv' - En + Es i{v + v' + lu) - En + El , 

X {ci)NMicl,)ML{Ca)Ls{Ca')sN (A. 3) 



q-PEm + e-PE_ 



^-PEs + ^-PE, 



E 

N,M,L,S 
o-PEa 



q-PEm _ ^-pEs 



q-PEl _ Q-PE^ 



i{v' + uj) + El- Em i{v + uj) + Es - En \ iuj + Es - Em i{v + v' + uj) ~ En + El 



-PE„ 



Xioc 



iv + Em - En iv' + El - Es 
-1 



X {cl,)NM{clr)ML{Ccr')Ls{Ca)sN 



(A.4) 



E 



N,M,L,S 



i{v' + Lu) + Em - E_ 



1 



N 



f^-pEL _ f^-pE 



^-PEm _^ q-0El q-0El _|_ q-PEs 



i{v - v') - Em + Es \iv + EL- Em iv' + El- Es 

X {cl,)NM{ci)ML{Ca')Ls{Ca)sN 



e~PE^ + ^-PEs 



i{v + uj) + Es - En \iiv + v' + uj) + El - En iv' + El - E, 



(A.5) 



Xioc 



E 



1 



i{v' +u;) + Em - En i{v + uj) + El - Es \ iuj + Em - Es iv + El - E_ 



N,M,L,S 
^-PEn _|_ q-PEs 



-PE^ 



, Q-PEs q~I3El _|_ q-PEm 



+ — 



-PEi, 



-PEl 



iv' + Es-En i{v + v' + Lu) + El - E. 



N 



X {cl,)NM{cl)ML{Ca)Ls{Ca')sN 



(A.6) 



Xioc 



E 



1 



N,M,L,S 



i{v' + uj)+Em - En 



1 



1 



-PEs _|_ Q-PE 



Q~PEs _ q~PEm 

i{v + uj) + Em - El \ iv + Es — El ' iuj + Em — Es 
e-PEs + ^-PE, 



i{v - v') + En - El \ iv + Es - El iv' + Es - En 

I 



X {cl,)NM(Ca)ML{cl)Ls{Ca')sN (A. 7) 



After the DMFT self-consistency condition has been fulfilled, the ED-DMFT evaluation of the local suscepti- 
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bility Eq. ((T7)) is obtained directly by plugging the eigen- 
values and the matrix elements {ca)NM of the asso- 
ciated AIM in the Eqs. (|A.2|l - (|A.7p and performing the 
corresponding four summations over the Hilbert space. 
The relevant numerical effort related to the Hilbert space 
summations can be handled by means of a parallel com- 
putations for the case considered here (i.e., TV^ = 5 and 
Nmax > 20). States for which all Boltzmann weights 
(e~^^) or matrix elements are smaller than 10"^ are ne- 
glected. 

For the numerical implementation, let us note that 



some denominators in Eqs. (jA.2p - (|A.7|) . characterized 
by a bosonic Matsubara frequency (e.g., those with iuj or 
with i{v — J/')), can vanish during the summations over 
the Hilbert space. 

However there are no divergences, since the corre- 
sponding numerators are simultaneously vanishing, so 
that their limit is always well defined. To avoid com- 
putational problems, we simply add a very small energy 
shift in all the terms with a vanishing denominator (e.g., 
i?jv — > ^'Af + 10~^), in order to evaluate numerically the 
correct limiting values. 
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